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EOARD  Grant:  FA8655-13-1-3052 
Report  Date:  March  26,  2015 

Final  report  for  the  duration  from  January  1,  2014  to  March  30,  2015 

1  Introduction 

The  Earth’s  ionosphere  is  the  well  known  threat  for  applications  that  depend 
on  trans-ionospheric  radio  waves.  For  example,  it  affects  applications  like  navi¬ 
gation,  positioning,  and  terrestrial  communication.  To  make  a  remedial  action 
for  the  effect  of  the  ionosphere  on  radio  wave,  the  behavior  of  the  ionosphere 
has  been  investigated  since  two-three  decades  ago  all  over  the  globe.  However, 
there  are  regions,  like  Africa,  that  have  not  been  well  investigated  until  recently 
due  to  lack  of  enough  ground-based  ionospheric  monitoring  devices.  In  addition, 
the  Low  Earth  Orbiting  (LEO)  satellites  ion  density  observations  show  unique 
features  for  the  African  sector  [Hei  et  al.  2005].  To  understand  the  unique  char¬ 
acteristics  of  the  African  ionosphere  different  organizations  within  and  outside 
Africa  are  being  participated.  For  example,  through  the  IHY-United  Nations 
Basic  Space  Sciences  Initiative  (UNBSSI)  a  number  of  ground-based  ionospheric 
monitoring  devices  have  been  installed  in  Africa  [Amory-Mazaudier,  et  al.  2009] 
since  2007.  Alongside  this  activity,  universities  in  Africa  (e.g.  Bahir  Dar  Uni¬ 
versity,  Ethiopia)  have  designed  space  science  curriculum  for  MSc  and  PhD  level 
and  they  are  training  young  scientists  in  collaboration  with  universities  and  re¬ 
search  centers  found  outside  Africa  (e.g.  ICTP,  Italy).  Washera  Geospace  and 
Radar  Science  Laboratory  (WaGRL)  at  Bahir  Dar  University,  is  a  good  example 
in  collaborating  with  other  institutions  (e.g.  AFRL,  DLR,  and  Kyoto  Univer¬ 
sity)  by  hosting  different  ionospheric  monitoring  instruments  and  carrying  out 
space  science  education  and  research  at  MSc  and  PhD  level. 

Similarly,  funding  agencies  are  supporting  the  activities  being  carried  out 
at  WaGRL.  The  EOARD  grant  FA8655-13-1-3052  is  the  best  example.  This 
final  report  summarizes  accomplishments  related  to  education  and  ionospheric 
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monitoring  facilities  and  research  findings  that  have  been  obtained  in  one  year 
research  contract  (  January  1  2014  to  March  30,  2015)  awarded  by  EOARD. 
In  the  coming  sections  we  describe  briefly  about  project  objectives  and  major 
achievements  and  research  findings  that  have  been  obtained  in  the  past  one  year 
and  three  months. 


2  Motivation  of  the  project 

LEO  satellite  observations,  in  the  African  sector,  show  unique  equatorial  iono¬ 
spheric  structure  [Hei  et  al.  2005].  For  example,  this  region  equatorial  plasma 
bubbles,  which  produce  scintillation  of  radio  waves,  are  much  deeper  and  occur 
more  frequently  than  bubbles  observed  in  the  other  region.  However,  these  ob¬ 
servations  have  not  been  confirmed,  validated  or  studied  in  detail  from  ground 
due  to  lack  of  ground-based  instruments  in  the  region,  and  the  question  what 
causes  these  unique  density  irregularities  in  the  region  is  not  yet  understood. 
Moreover,  global  empirical  models,  which  are  useful  for  ionospheric  study  and 
related  applications,  have  been  developed  using  ionospheric  data  mainly  col¬ 
lected  from  the  Northern  mid-latitude  regions.  This  implies  that  these  models 
are  producing  the  African  ionospheric  characteristics  by  interpolation  and  it  is 
worthy  to  investigate  answers  for  the  following  question:  how  do  the  existing 
models  perform  in  estimating  and  forecasting  the  unique  ionospheric  structures 
in  the  African  sector?  To  get  answers  for  these  and  similar  questions  (in  general 
to  understand  and  model  the  African  region  ionosphere)  space  science  educa¬ 
tion  and  research  are  being  carried  out  at  MSc  and  PhD  level  at  Bahir  Dar 
University.  In  addition  to  education  efforts,  new  coherent  scatter  radar,  which 
will  be  the  main  data  source  for  equatorial  irregularity  study,  has  been  installed 
in  collaboration  between  AFRL  and  Bahir  Dar  university.  These  activities  that 
have  been  carried  out  at  Bahir  Dar  University  (mainly  by  PhD  students)  are 
useful  to  complement  the  efforts  that  are  being  carried  out  ,  in  particular,  to  un¬ 
derstand  the  unique  features  of  the  African  equatorial  ionosphere  and  in  general 
to  understand  the  behavior  of  the  ionosphere. 

3  Project  achievements  and  research  findings 

3.1  Coherent  Backscatter  Radar  installation 

WaGRL  houses  different  ionospheric  monitoring  instruments  such  as  SCINDA 
GPS  receiver,  Galelio  receiver,  LEO  satellite  signal  receivers,  VLF,  VHF  and 
so  on  in  collaborating  with  institutions  found  outside  Ethiopia.  As  a  continua¬ 
tion  of  hosting  instruments,  HF  Coherent  Backscatter  Radar  named  Blue  Nile 
Radar  is  installed  at  Bahir  Dar  University  in  collaboration  with  AFRL  during 
November  2014.  The  radar  antenna  and  sample  result  are  shown  in  Figure  1. 
Installation  of  the  radar  is  one  of  achievements  of  the  project  supported  by 
EOARD  Grant  FA8655-13-1-3052.  Data  from  this  system  will  be  analyzed  for 
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further  study  to  identity  the  cause  for  unique  feature  of  the  African  equatorial 
ionosphere. 


a 


Figure  1:  Coherent  Backscatter  radar  antenna  (a)  and  SNR  as  a  function  of 
range  and  time  (b)  obtained  on  14  February  2015. 


3.2  Education 

Space  science  education  is  being  carried  out  at  PhD  level  at  Bahir  Dar  university 
to  produce  experts  who  will  contribute  for  the  understanding  of  space  physics. 
Supported  by  EOARD  grant  FA8655-13-1-3052  six  PhD  students  have  been 
studying  space  science.  In  the  history  of  physics  department  at  Bahir  Dar 
university  the  first  PhD  student  (Melessew  Nigussie)  in  space  science  (supported 
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Figure  2:  PhD  graduation  photo. 


by  TWAS,  ICTP  in  Italy  and  EOARD  grant  FA8655-13-1-3052)  has  obtained 
his  PhD  on  July  2,  2014  (see  graduation  photo  in  Figure  2).  His  PhD  work  has 
been  about  modeling  the  East  African  ionosphere.  Dr  Melessew  has  graduated 
after  publishing  three  papers  in  well  known  journals  (e.g.  AGU,  Radio  Science). 

3.3  List  of  publications 

During  this  project  time,  we  have  made  significant  research  findings,  which 
address  some  of  the  research  questions  listed  above.  The  results  are  found  in 
the  papers  and  manuscripts  listed  below. 

•  Nigussie,  M.,  B.  Damtie,  E.  Yizengaw,  S.  M.  Radicella,  Modeling  the  East- 
African  ionosphere,  America  Geophysical  Union  (AGU)  Book  Chapter,  in 
press,  2015. 

•  Asmare,  Y.,  Kassa,  T.,  Nigussie,  M.,  Validation  of  IRI-2012  TEC  model 
over  Ethiopia  during  solar  minimum  (2009)  and  maximum  (2013)  phases, 
Advances  in  Space  Research,  2014- 

•  Kassa  T.,  Damtie  B.,  Bires  A.,  Yizengaw  E.,  Cilliers  P.,  Storm-time  char¬ 
acteristics  of  the  equatorial  ionization  anomaly  in  the  East  African  sector, 
Adv.  Space  Res.,  on  review. 

•  Kassa  T.,  Damtie  B.,  Bires  A.,  Yizengaw  E.,  Cilliers  P.,  Spatio-temporal 
characteristics  of  the  equatorial  Ionization  Anomaly  (EIA)  in  the  East 
African  region  via  ionospheric  tomography  during  the  year  2012,  Advances 
in  space  research,  55,  184-198,  2015. 

•  Tesema,  F.,  B.  Damtie,  M.  Nigussie,  The  response  of  the  ionosphere  to 
intense  geomagnetic  storms  in  2012  using  GPS-TEC  data  from  East- Africa 
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longitudinal  sector,  submitted  for  publication  to  Atmospheric  and  Solar- 
Terrestrial  Physics  journal ,  it  is  on  review. 

•  Tebabal  A.,  B.  Damtie,  A.  Bires,  M.  Nigussie,  E.  Yizengaw,  Reconstruc¬ 
tion  of  total  solar  irradiance  using  artificial  neural  networks,  submitted  for 
publication  to  Atmospheric  and  Solar- Terrestrial  Physics  journal,  it  is  on 
review. 

3.4  Modeling  and  characterization  of  the  East- African  iono¬ 
sphere 

Ionospheric  empirical  models  such  as  NeQuick  and  IRI  have  been  developed 
using  data  collected  from  ionosonde  installed  unevenly  on  the  Earth.  Most 
of  the  data  used  for  these  model  development  have  been  collected  from  the 
mid-latitude  of  the  Northern  hemisphere.  This  implies  that  these  models  are 
being  used  to  estimate  ionospheric  characteristics  for  African  region  by  inter¬ 
polation.  So,  for  realistic  ionospheric  characteristics  estimation  for  the  African 
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Figure  3:  Diurnal  variations  of  the  contour  plot  of  electron  density  profile  esti¬ 
mated  using  NeQuick  driven  by  Az  maps  at  39°  E  on  March  11,  2011 


ionosphere  these  models  should  be  assisted  by  measurements  collected  recently. 
For  near-real  time  ionospheric  characteristics  (electron  density  and  total  elec¬ 
tron  content  (TEC))  estimation,  new  data  driven  modeling  technique,  which  is 
used  to  combine  the  GPS  TEC  and  the  existing  empirical  model  (NeQuick), 
has  been  developed.  Sample  results  from  this  method  are  shown  in  Figures  3, 
4,  and  5.  Figure  3  shows  the  diurnal  variations  of  the  vertical  electron  density 
distribution  estimated  at  39°E  meridian.  The  plots  show  clear  spatiotemporal 
variation  of  electron  density.  Evidently,  in  the  pre-and  post-sunset  hours  (1600 
to  2100  UT)  symmetric  and  clear  EIA  anomalies  are  depicted  unlike  that  of 
the  results  obtained  between  1000  and  1500  UT.  This  indicates  that  the  plasma 
fountain  effect  in  between  1600  to  2100  UT  is  much  stronger  than  the  plasma 


5 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


fountain  effect  between  1000  to  1500  UT,  which  could  be  due  to  strong  pre¬ 
reversal  enhancement  of  ExB  [Fejer  et  al.  1991]. 

The  methods  developed  in  this  work  is  also  useful  to  compute  three  dimen¬ 
sional  electron  density  profiles.  Figure  4  shows  the  three-dimensional  electron 
density  maps  estimated  at  1700  UT.  This  plot  shows  the  latitudinal,  longi¬ 
tudinal,  and  altitude  variations  of  the  electron  density  estimated  after  model 
adaptation.  As  can  be  seen  in  this  figure  while  clear  latitudinal  and  altitude 
variability  are  visible  there  is  no  that  much  significant  longitudinal  variation  in 
the  electron  density  distribution,  indicating  the  ionospheric  density  vary  less  as 
a  function  of  longitude. 


Three  dimesional  electron  density  distribution  at  1700  UT 


Altitude  [Km]  Geog.  La,  [Degree] 


Figure  4:  Three  dimensional  electron  density  maps  estimated  by  non- 
tomographic  technique  at  1700  UT  on  March  11,  2011. 


In  order  to  validate  the  effectiveness  of  the  non-tomographically  produced 
density  maps,  conjunction  C/NOFS  trajectories  have  been  identified.  Fig¬ 
ure  5  shows  the  trajectories  of  the  C/NOFS  satellite  (top  panels)  at  different 
local  time.  The  corresponding  experimental  (from  PLP  instrument  onboard 
C/NOFS)  and  modeled  (before  (F10.7)  and  after  (Az)  data  adaptation)  elec¬ 
tron  densities  are  shown  in  the  bottom  panels  in  Figure  5.  As  expected  the 
in-situ  density  estimated  after  the  model  adaptation  agree  very  well  with  the 
C/NOFS  observations,  except  when  the  satellite  was  around  F2  peak  altitudes 
(  400km)  as  shown  in  the  right  bottom  panel  in  Figure  5.  This  discrepancy 
might  be  due  to  the  slight  time  difference  of  the  ionospheric  information  uti¬ 
lized  for  model  adaptation  (at  1000  UT)  and  obtained  by  the  C/NOFS  satellite 
(between  10.0353  and  10.1254  UT).  The  method  that  has  been  used  to  produce 
these  results  have  been  published  recently  in  AGU  books  [Nigussie  et  al.  2015]. 
This  method  will  also  be  applied  to  investigate  different  ionospheric  phenomena 
such  as  the  effect  of  the  geomagnetic  disturbance  on  the  distribution  of  the  low 
latitude  ionospheric  electron  density. 

Similarly,  ionospheric  tomography  method  has  been  applied,  starting  from 
basic  principles,  to  reconstruct  the  ionospheric  electron  density  from  the  GPS  re- 
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C/NOFS  trajectory  (4.8543  to  4.9476  UT) 


C/NOFS  trajectory  (10.0353  to  10.1254  UT) 


550 


x  1011 
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Geo.  Lat.  [Degree]  Geo.  Long.  [Degree]  Geo.  Lat.  [Degree]  Geo.  Long.  [Degree] 


Figure  5:  Trajectories  of  C/NOFS  satellite  (top  panels)  and  modeled  (marked 
by  Az  and  F10.7)  and  experimental  (C/NOFS)  electron  densities. 


ceiver  measurements.  The  MATLAB  code  developed  to  implement  this  method 
is  one  of  the  outputs  of  this  project.  Sample  reconstructed  electron  density 
maps  using  this  method  are  shown  in  Figure  6.  The  results  show  the  spatio- 
temporal  evolution  of  the  electron  density  reconstructed  from  the  GPS  receiver 
measurements  taken  in  the  East- African  ionosphere.  This  result  and  other  out¬ 
puts  of  the  method  are  published  [Kassa  et  al.  2015].  Analysis  of  data  has  also 
been  carried  out  to  study  the  effect  of  the  geomagnetic  storm  on  the  equatorial 
ionosphere.  Sample  result  is  shown  in  Figure  7.  This  figure  shows  the  spatio- 
temporal  variation  of  TEC  (top  two  panels)  and  equatorial  electroject  (bottom 
panels)  strength  (in  nT)  estimated  from  magnetometers  on  07-11  March  2012. 
From  the  study  that  has  been  carried  out  to  get  this  plot,  we  have  found  out 
that  the  variation  of  TEC  for  East- African  longitudinal  sector  mainly  affected 
by  the  modification  of  low  latitude  electrodynamics  and  neutral  composition 
change  due  to  the  storm  effect.  The  work  that  has  been  carried  out  to  get  this 
sample  result  has  been  submitted  for  publication  to  Atmospheric  and  Solar- 
Terrestrial  Physics  journal  and  now  it  is  on  review. 

The  equatorial  ionosphere  is  known  on  its  effect  on  radio  waves.  Rapid  fluc¬ 
tuation  of  the  ionosphere  severely  affects  a  GPS  signal.  The  fluctuation  can  be 
described  by  the  scintillation  index,  which  can  be  estimated  from  received  power 
of  a  GPS  satellite  signal  or  from  the  rate  of  change  of  total  electron  content. 
This  means  that  it  is  important  to  have  good  approach  to  estimate  scintillation 
index  (S4)  from  the  TEC  measurements.  During  this  project  period,  stochas¬ 
tic  approach  has  been  applied  to  estimate  the  rate  of  change  of  TEC  (ROTI). 
Sample  output  of  the  method  applied  is  shown  in  Figure  8.  The  figure  shows 
the  TEC  variation,  S4  index  estimated  from  SCIND  data,  and  ROTI  estimated 
from  TEC  using  the  stochastic  approach,  from  top  to  bottom,  respectively.  As 
can  be  seen  from  this  result,  ROTI  from  stochastic  modeling  and  S4  index 
from  SCIND  A  system  equivalently  describe  the  fluctuation  of  the  TEC.  This 
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Figure  6:  Temporal  variations  of  two  dimensional  electron  density  distribution 
for  the  East- African  ionosphere. 


approach  is  still  under  development  and  hence  a  manuscript  will  be  sent  to  a 
journal  for  publication. 

In  addition  to  ionospheric  study,  we  have  been  involved  to  model  the  characteristics 
of  the  solar  variability,  which  effects  directly  the  variability  of  the  ionosphere.  In 
this  modeling  technique  we  have  applied  the  neural  networks  package  available 
in  MATLAB  toolbox,  to  represent  the  trend  of  the  data  and  use  it  for  prediction. 
Sample  results  are  shown  in  Figure  9.  The  results  show  the  comparison  of  the 
total  solar  irradiance  (TSI)  output  from  neural  networks  after  data  assimilation 
(green  dots)  with  that  of  the  measurements  (black  line)  that  do  not  involved  in 
assimilation.  The  comparison  shows  excellent  agreement.  This  and  similar  re¬ 
sults  have  been  submitted  for  publication  in  Atmospheric  and  Solar- Terrestrial 
Physics  journal  and  now  it  is  on  review. 
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Figure  8:  Variation  of  TEC  (top  panels),  S4  (middle  panels),  and  ROTI  (bottom 
panels). 


Figure  9:  Reconstructed  (green)  and  measured  (black  line)  TSI. 


10 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


